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Abstract 

The longitudinal structure function of the (i(e, e'p) exclusive cross section is 
calculated with the Lorentz integral transform method. In this approach final 
state interaction is fully taken into account, but without using a final state 
wave function. Cross sections are obtained via the inversion of the transform. 
It is shown that the inversion results are very stable. The comparison to a 
conventional calculation with an explicit np final state wave function shows 
that the obtained results are also very precise. Thus the method opens up 
the possibility to obtain exclusive cross sections for reactions with more than 
two particles, where it is generally very difficult to calculate the exact final 
state wave function. 
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I. INTRODUCTION 



The study of inclusive and exclusive cross sections in inelastic reactions is an essential 
tool in understanding the underlying dynamics of a particle system. For systems with more 
than two constituents a major problem in the calculation of such reactions consists in the 
exact knowledge of the final state wave function in the continuum. Its calculation is far 
more difficult than the corresponding bound state calculation. As a matter of fact, today a 
computation of an intermediate energy continuum wave function is out of reach for a system 
with more than three particles. An exact calculation, however, can be carried out in an 
alternative way by using proper integral transforms. They allow one to take into account 
final state interactions (FSI) rigorously without using final state wave functions explicitly. 
In fact in recent years the Lorentz integral transform (LIT) method [|I| has been successfully 
applied to various inclusive breakup cross sections in few-body physics. The LIT is expressed 
in terms of square integrable functions which are obtained from inhomogeneous differential 
equations. The differential equations can be solved with similar methods as a bound state 
problem. After having calculated the transform one obtains the cross section by inversion 
of the transform. A recent overview of the results obtained with the LIT method is given in 
i- 

The LIT method has not yet been used for the calculation of exclusive cross sections. 
The only exception is a calculation of the ^He spectral function [Q, but in this case the 
method proceeds along the same lines as for an inclusive reaction. However, it is in principle 
possible to apply the integral transform method to general exclusive reactions as shown by 
Efros for the Stieltjes transform |Q. In the present paper we illustrate the details of the 
calculation for the d{e', ep) reaction with the LIT. It will also allow us to check the precision 
of the obtained results by comparing them to those of a conventional calculation with an 
explicit np final state wave function. Such a reliability check is necessary, since the use 
of integral transforms is not always unproblematic for the calculation of cross sections. In 
fact one may encounter problems in the inversion of the transform. For example in Ref. 
it was shown for the case of an inclusive reaction that the Stieltjes transform is not very 
appropriated, since it samples contributions over a large energy range making the inversion 
extremely difficult. Later the LIT was proposed and it was shown for the test case of the 
longitudinal inelastic deuteron form factor in electron scattering that inclusive reactions can 
be safely calculated with this method. Whether the LIT is as appropriated also for exclusive 
reactions cannot be said a priori, since such a calculation is more complicated. The aim of 
the present work is to investigate this question. 

The paper is organized as follows. In Sec. || we describe how cross sections are calcu- 
lated with the LIT. Details of the calculation for the specific reaction under consideration 



{d{e, e'p)n) are given in Sec. |T|. The results are illustrated in Sec. and a conclusion is 
drawn in Sec. |V|. 

II. THE LIT METHOD 

The starting point of the LIT method [ij is the calculation of an integral transform with 
a Lorentz kernel 
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L{a) = I dco- 2 • (1) 

The function F depends on the internal excitation energy u = Ef — Eq of a given particle 
system and contains information about the transition of the system from the ground state 
I^E'o)) with energy Eq, to the final state l"^ /), with energy Ef, induced by an external probe. 
In case of an inclusive reaction F{uj) denotes the response function 

F{uj) = J d<i/f\{^f\d\-^o)\'^S{Ef - Eo-uj) , (2) 

where O is a transition operator which characterizes the specific process under consideration. 

The key point of the LIT method is an evaluation of L{a) without explicit knowledge 
of F{uj). In a second step the function F is obtained from the inversion of the transform. 
The great advantage of the method lies in the fact that a calculation of the generally very 
complicated final state wave function |\E'/) can be avoided as will be discussed below. On 
the contrary a conventional calculation of F{uj) can only be carried out with the explicit 
knowledge of I^P/). 

Using completeness one can show that in order to obtain L{a) one has to solve the 
following differential equation 

{H-Eo~a*)\^ii(T)) = d\^o) (3) 

with 



(J = ctr + iai aR, aj >0, (4) 

where H = T + V is the Hamiltonian of the system under consideration. Note that the cor- 
responding homogeneous equation has only the trivial solution, since H has a real eigenvalue 
spectrum. 

The norm of the solution \E'i(cr) determines the LIT directly: 

L(a) = (^i(a)|^i(o-)). (5) 

Different from a Schrodinger equation at positive energies, one has for the solution of Eq. 
(0) a very simple boundary condition. Due to the localized source at the right hand side 
(rhs) of Eq. @ ^i{<y) vanishes at large distances similar to a bound state wave function. 
Therefore one can apply the same techniques as for the calculation of a bound state wave 
function. 

The response function F{uj) serves only for the determination of inclusive cross sections. 
For an exclusive process one needs a more detailed information about the transition of the 
system. In fact one has to be able to calculate transition matrix elements of the form 

Tf,{Ef) = {^f\d\^o)- (6) 

How such a calculation can be carried out with an integral transform was shown for the case 
of the Stieltjes transform For the LIT the calculation proceeds analogously as outlined 
in the following (see also ||,^). 
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For simplicity we will consider an exclusive reaction leading to a final state with two 
fragments, but the method can be applied also for channels with more than two particles. 
Besides the correct final state wave function we also introduce the corresponding plane 
wave 

4™(f) = A^.'i/^^^g^ . (7) 

where ^ is a proper antisymmetrizer and and \l/2 are the internal wave functions of the 
two fragments, while k and f are the usual relative coordinates for momentum and position 
of the two-body system formed by the two fragments. 

In order to calculate the T matrix element one starts from the Lippmann-Schwinger 
equation for the final state 

with 

V = J2V,,, zGJ-i,jG^2, (9) 

where J^i and J^2 contain all particles of the first and second fragment, respectively. Inserting 
the above expression in Eq. (|^) one obtains a sum of two pieces, a Born term 

Tfr{Ef) = {'i>^'^\d\^o), (10) 

and a term depending on FSI 

Tj^Ef) = (<f^^|V^--^ -d\^o) ■ (11) 

tjf + ie — rl 

The calculation of the Born term is rather simple. The main difficulty of the calculation is 
the determination of the matrix element depending on FSI. 

Using the completeness of the eigenstates \^{E)) of H one can rewrite Tji'^ as follows 

Tj!\E,) = |dB{4™|('|*(B))(*(B)|_l--^6|*,) 

^E^d^)!^;^^ . (12) 

with 

Ff,{E) = {^'''^\V\m{E)){m{E)\d\^>,) . (13) 
One obtains the following formal solution for the FSI term 

Tf\Ej) = -^nF,,{E,) + V Q dE§^ • (14) 
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For a calculation of Tji one needs to know Ffi for any given energy. A direct calculation of Ffi 
is of course in general far too difficult, since one has to determine final state wave functions 
\l/ for the whole eigenvalue spectrum of H. On the other hand an indirect calculation via 
the LIT is possible. To this end one performs a Lorentz integral transform of Ffi, i.e. 

Inserting the definition of Ffi one finds 

L(a)= / dE{<l>^'^\V- mE)){nE)\j, -0\^o) 

JEq H — a H — a 

H — a H — a* 

= {^,{a)\^,{a)), (16) 

with 

l^i(^))= ^ \- 0\^o), (17) 
\^2{cy)) = j^-^——V\^^'^). (18) 

It is evident that (|l^ leads essentially to the same differential equation as for the inclusive 
process: 

{H-an + iai)\^r{a)) = d\^,). (19) 

From Eq. (|18]) one obtains 

{H-aR + iai)\i>2{a)) = f l^^^) , (20) 



which is similar to Eq. (IS), but for a different source term on the rhs. For a finite 
range potential one has also in this case a vanishing source term for large distances. This 
guarantees also for an asymptotic boundary condition similar to a ground state problem. 
If V contains also the Coulomb potential one cannot proceed exactly in the same way as 
shown here. In this case one has to start from a modified Lippmann-Schwinger equation, 
where Coulomb wave function are taken into account 

As shown by Efros j^,]^ the integral transform method can be extended to exclusive 
processes with more than two fragments. In principle one obtains equations similar to Eqs. 
(^,|T^,^). However, one cannot guarantee, as in the two fragment case, that the potential 



V vanishes asymptotically, since two of the fragments could remain close to each other. 
Therefore it is necessary to choose a different solution to the problem. In fact one can 
rewrite L{a) as follows 

L(a) = (<l>PW|F|l'i), (21) 
where |\E'i) is obtained from the solution of the following differential equation 
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[H-aR- iaj)\^iia)) = |^i((t)) . (22) 

Since |\E'i) vanishes at large distances, one can again use bound state methods for the solution 
of this differential equation. 

There is another possibility to calculate the LIT M . Starting from the following identity 
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1/1 In 

^r\H - n ~ H - rr*) ' 



and defining 



{H-a){H-a*) 2iai^H - a H - a* 



m^)) = j7 ^^Ol^o), (24) 



one gets for the transform 

L{a) = ^^{{^""^iVl^'ii^)) - ($^^|V^|^i(«-))) . (25) 

It is seen that one has to solve only one type of differential equation. In fact for the solution 
of |\^'^(cr)) it is sufficient to solve Eq. (|19D a second time replacing a* by a. 

Which of the two approaches for a nonvanishing source term is more appropriate depends 
also on the possibility to obtain a precise numerical solution. Small errors in the solution of 
Eq. ( p!9D could lead in both cases to larger errors for the calculation of the LIT. Therefore 



we will also consider in Sec. [V these additional possibilities to calculate the LIT for our 



realistic test case of the electromagnetic deuteron breakup. 



III. EXCLUSIVE DEUTERON BREAKUP 

The exclusive deuteron breakup d{e, e'p)n is governed by the four structure functions 
fr, fiT, and fxr- In the following we will only consider the longitudinal Z^. We use the 
same notation as in Ref. |0, where the structure functions are calculated in the final np cm. 
system with fi = fii^Enp, \qc.m.\,0)- In this system qcm. denotes the momentum transfer, 
Enp the relative np energy, and the relative np momentum k has the angle 6 with respect to 
qcm.- The structure functions are expressed in terms of the transition matrix T^miimi for the 
process e + d ^ e' + np. The quantum numbers S and m denote the spin and spin projection 
of the outgoing np pair with respect to k, /i characterizes the transition operator O, and 
is the projection of the deuteron spin with respect to qcm.- For the longitudinal structure 
function, i.e. /i = 0, one has the following transition operator 

O = X! ^eMI) exp(ig;.m. ■ r^) , (26) 
j 

where Ge,] denotes the electric form factor of the j-th nucleon with being the four- 
momentum transfer squared. Here we use the electric dipole form factor for the proton, 
while the neutron electric form factor is set to zero. It is evident that the above operator 
does not affect the spin, i.e. 5* = 1, and one gets ^ 
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where 

TimOma = C(lm|6|mrf) (28) 

with 




C = -{2n)^-—, (29) 



where M denotes the nucleon mass and a is the fine structure constant. Like fi the transition 
matrix depends on E^p, qc.m., and 9. In the following we suppress the dependences on gc.m. 
and 6, while the dependence on Enp is made explicit in most cases. 

The aim of the present work is a test of the LIT method for the exclusive deuteron 
breakup. Since there is no need to obtain results with a realistic potential we choose the 
semi-realistic TN potential for this test. It is a central potential model, which is different 
for the various spin-isospin (ST) channels, i. e. V{r) = Z^st ^^"'"('")- Since we have 5 = 1, 
only V^^ and V^^ have to be considered here. In addition V^^{r) = for the TN potential, 
hence FSI effects appear only in the channel ST=10. 

Due to the absence of the tensor force in the TN potential one obtains a reduced com- 
plexity of the equations to be solved. Nevertheless, as calculations with this potential model 
show (see e.g. 0) results are sufficiently realistic in order to serve as test case. 

As discussed in Sec. |T| there are two contributions for any T-matrix element, a Born and 
an FSI term. Considering an s-wave deuteron with radial wave function u{r) one obtains 
for the Born term 

Ti^om, = i-r''6^^,GE,pV3kMYooik_) H drruir)Uk_r) , (30) 

Jo 

where jo denotes the spherical Bessel function of order and 

k^ = k-^^. (31) 

For the second piece, T^^\ it is necessary to perform a multipole decomposition. It is 
convenient to introduce the following expansions with projections M with respect to qc.m. 

^iM = Y.^'^[Y^'^Kr)xx^'\a,,a2)filV2LTl (32) 
^2,M = E ^'[Y^Hr) X xW(^i, ^2Mcl^,mMYLfi)r-'^l-\r) , (33) 

j,l,mi 

where x'^^(^i; ^2) denotes the spin wave function for a two-nucleon system with 5* = 1. For 
the rhs of the differential equations (0) and ( pOD we perform similar expansions as for 
and leading to 

0|m,) = E^'^[>^f^kOxx''k^i,^2)]l^l,v/2LTTC'o^J.^„,r-VL(r), (34) 
hir) = jLi^-^Mr) (35) 
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and 



V^I^D = E X X^H^i^^2)]W^XMyLfiy''9iAr) , (36) 

j,l,mi 

9ijir) = ]j^rji{kr)Vji{r) , (37) 

where Vji is the potential for the NN partial wave ^Ij. 

For the above multipole decompositions lead to the following coupled differential 
equations in real and imaginary parts 

{-^ - + V^r) - a«}5R[^f;V,r)] - aj-s[^pif {a,r)] = g,,{r) (38) 

{-| fl^ - ^1 + ^'^-^^^ - + r)] = . (39) 



As mentioned above, because of our potential model we have to consider only the channel 
with ST=10, thus Vij{r) can be replaced by V^^{r). In addition the Pauli principle has to be 
fulfilled, i.e. S+T+/ has to be odd. Therefore the differential equation has only to be solved 
for I even. For channels with 1 odd one had to consider V^^, but as already mentioned V^^ 
is zero in our potential model. Note that there is no explicit dependence on j in the coupled 
differential equation, thus one has 'iplj^ = ipij) . 

For \l'i one finds very similar equations with the only difference that one has to replace 
gij{r) by fL{r) on the rhs. Also here we have ip^j^j = ipij'- 

We solve the differential equation by adding an additional homogeneous equation deter- 
mining the source terms on the rhs. In this way we obtain a coupled homogeneous differential 
equation system. The numerical solution leads to very precise results as shown in Ref. [^. 

With the solutions for ijj^j^j and ip^j^ one obtains for the scalar product ([T6|) 

(^2|^l)= E V2rTlYUk)T.Cl^UrnA'^.rnM^) (40) 
mnialmi j 

with 

Ri,{a) = dr{^lf{a, r))>«(a, r)5u . (41) 
To calculate the FSI contribution to the T-matrix elements one has to invert the LIT 

in order to obtain the function Rij{E). The transform can be inverted using the following 
ansatz 

N 

RijiE) = J2cn,i,Xn,ijiE,P), (43) 

n=l 

where Xn,ij are given functions with nonlinear parameters jS. Substituting this expansion 
into the rhs of Eq. (^21) one obtains 
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N 

RlM) = Cn,ljXn,lji(^, P) , (44) 
n=l 

where the Xn,zi are the Lorentz integral transforms of the basis functions. The parameters 
Cn^ij and (3 are determined by fitting the calculated transform Rij{a) to the above expansion 
at many points for a fixed aj. The number of functions plays the role of a regularization 
parameter and is chosen within a stability region, i.e. where the obtained results are stable 
for a certain range of N (see also 0,0]). Here we use the following set of basis functions 



X^,,iE,/3) = E'+hxpi-^). (45) 



For the parametrization of the elastic monopole transition for / = we include an additional 
function in the set x'- 

Xo,oAE,P) = 6{E-Eo). (46) 



Thus the sum in Eqs. (^, ^4[) starts in this case with n = instead of n = 1. 

Once the inversion is carried out one can make use of Eq. (|T^). For a specific set of m 
and rrid one finds the following FSI contribution to the T-matrix 



TlmOma i^np) — X! (^) l/ ^ImOmd i^np) (47) 



Imi 



with 

For our case without tensor force the following simple relations hold 

,FSI,« _ r ,FSI,« 

,FSI,« _ .FSU _ .FSl,l fro\ 
t-l-lO-l ~ ''1000 ~ ''llOl • K^'^) 

It is worth mentioning that as a byproduct of the calculation one obtains also the NN 
phase shifts from our calculation, i.e. without having solved the Schrodinger equation for the 
scattering state. The ratio of imaginary and real parts of a given transition matrix element 
is equal to tan{S) (see e.g., Ref. 0). For our simple potential model one obtains 

5i = atan( ^ , ,, 5 — r ), (51) 

where is analogously defined as t^^^'^, and easily evaluated from a multipole decompo- 

sition of the Born term in Eq. (0). If one is only interested in the phase shifts themselves, 
one can perform a simpler calculation neglecting the excitation operator O (see Ref. M). 
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IV. DISCUSSION OF RESULTS 



We test the LIT method for exclusive reactions choosing for the electromagnetic deuteron 
breakup three different kinematics with rather strong FSI effects: (i) in the tail region beyond 
the quasi-elastic peak at moderate momentum transfer {Enp = 120 MeV, = 5 fm~^), (ii) 
on the photon line {E^ = 70 MeV), and (iii) close to the deuteron breakup threshold {Enp = 1 
MeV, Qcm ~ fm~^). Note for kinematics (ii) that there is no longitudinal contribution to 
the exclusive (e, e'p) cross section, but that nevertheless the structure function does not 
vanish. In fact in applying Siegert's theorem the longitudinal matrix elements are commonly 
used in photodisintegration and lead to the dominant contribution for the electric transitions. 

In order to have a more detailed comparison between the LIT results and the results 
of a conventional calculation we do not simply discuss the final result for fi, but rather 
study directly the FSI effect on the various multipole transitions. This allows us to make a 
much more precise comparison between the two calculations. Because of Eqs. (EDLEDI) it is 

FSI / I— r 

sufficient to consider tuQi in the following. 

Before turning to the above mentioned three kinematical cases we first illustrate results 
for the transform in a more general way choosing a constant of 5 fm~^ and various 
energies Enp. In Fig. 1 we show Rii{aji,ai = 20 MeV) for Z = 0, 2. Its inversion, Ru, 
gives a contributes to t^si,/ ^^^^ (^5])). For the / = transforms one has an interesting 
structure at small aji. It originates from the rather strong monopole transition strength 
close to the deuteron breakup threshold. It is interesting to see that the peak becomes more 
and more pronounced for the case that also Enp moves closer to the threshold region. In 
addition there is a second rather sizable contribution in the low ctr range. It arises from the 
elastic monopole contribution. Therefore one has to pay attention in the inversion of Rqq. 
One has to check whether aj is small enough to resolve with sufficient precision the elastic 
contribution from the threshold contribution. For I = 2 one has a rather different picture. 
One finds a peak in the quasi-elastic region. Note that for the considered momentum of 
q'^ni. of 5 fm~^ the quasi-elastic peak is situated at about Enp = 50 MeV. To find such a 
quasi-elastic peak for R22 is a bit surprising, since one does not expect there strong FSI 
effects. On the other hand FSI should be small not on an absolute scale but compared with 
the corresponding Born term. Furthermore, the real part of the FSI contribution is difficult 
to estimate from Fig. 1 because a principle value integral has to be calculated in this case 
(see Eq. 

In Fig. 2 we show Rii{E) of kinematics (i) for Z = 0, 2, 4 and cr/ = 5 and 20 MeV. It is 
obtained from the inversion of the corresponding Ru^aR^aj) (see Eq. (^2])). One sees that 
six basis functions are not sufficient for the inversion, but for a higher one obtains a very 
nice stability of the inversion. Comparing the results with different o"/, one finds a small 
difference for / = in the threshold region. The differences arise because the monopole 
contribution has a peak at the very threshold which has to be separated from the elastic 
contribution at Eq = —2.225 MeV. From the inversion we obtain an elastic contribution of 

3 

about 1.1 fma which is rather sizable compared to the inelastic part with a peak height of 

3 

0.045 fm2. Therefore is not surprising that the higher resolution with o"/ = 5 MeV leads 
to a somewhat different result. However, because of the rather high Enp of 120 MeV, the 
difference at the threshold is rather unimportant. This is confirmed by the results for t^^^, 
which are shown in Fig. 3 as function of the number of basis functions used for the inversion. 
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In fact the agreement among the results with o"/ = 5 and 20 MeV is very good. It is seen 
that one obtains for > 10 for all considered multipolarities I and for both ct/ values very 
similar and stable results. 

Also shown in Fig. 3 is the t^^^ of a conventional calculation. These results are very 
similar to the LIT results with relative differences of less than 1%. Only for the real part 
of the I = 4 transition the difference is a little bit larger. On the other hand one has 
also to consider that this matrix element is very small. In fact its size is only —2.5% of 
the corresponding t^'^'^^ matrix element. Thus the relative difference for the total matrix 
element is of the order of 10~^. For such a small FSI effect a part of the differences could 
also be due to a not completely exact result of the conventional calculation. Different from 
the LIT method t^^^ is not calculated directly, but taken indirectly from the difference of 
^totai g^j^j t^°"^-^ here tt°*^i corresponds to the transition with the correct np final state wave 
function in presence of the potential. The FSI effect is much more sizable for the two other 
transitions. Taking also here the ratio of ^Fsy^Sorn f^j. ^j^g j-g^^^ parts, one finds for Z = 
about -40% and for / = 2 about -30%. 

Results for kinematics (ii) are shown in Figs. 4 and 5. For the Ru{E) of Fig. 4 one 
finds again nice stabilities of the inversion for a larger number of basis functions. Comparing 
the Rii with the two different o"/ one has also here differences for the monopole transition 
and in addition for I — A. The monopole is of course not relevant for this kinematics on 
the photon line, since there is no corresponding electric monopole. On the other hand it 
is interesting to see whether one is able to separate the strongly peaked threshold strength 
from the dominant elastic contribution. Due to the lower momentum transfer one obtains 
an even larger elastic Ru than for kinematics (i), namely a value of about 14 fm^ . The 
^FSi j-gg^h;s QT^Q shown in Fig. 5. The real parts turn out to be very stable as function of 
number of inversion basis functions. They are also very similar for both o"/. Here we have 
the following relative FSI effects comparing with the Born term: —100% {I — 0), —25% 
(1 = 2), -1.5% (/ = 4). 

In comparison to the conventional calculation one finds in Fig. 5 for all the real parts 
only very small differences of less than 1%. For the imaginary part of t^^^ the picture is 
a bit different. The I = 2 results are very stable and agree with extremely high precision 
to the results of the conventional calculation. Also the / = 4 results are stable, but they 
are a few percent larger than found in the conventional calculation. However, one should 
note that the matrix clement is very small and hence the difference of a few percent is not 
relevant. In fact comparing with the above mentioned size of the real part of the total matrix 
element the difference between both calculations is of the order of 10~^. The imaginary part 
of the I = transition shows a bit less stability with the number of inversion basis function 
rcficcting also the above mentioned problems for the separation of the elastic contribution. 
On the other hand one obtains reliable results for the highest A^'s. 

For the third kinematics we illustrate the results in Figs. 6 and 7. Here we consider only 
I = and I = 2 transitions, since FSI effects do not play any role for higher transitions at 
threshold. In fact the FSI contribution is already very small for Z = 2. The inversion results 
in Fig. 6 are again very stable, except for / = with o"/ = 20 MeV. Of course, again it is the 
problem associated with the elastic contribution (i?oo(-^o) is about 15 fm^). Figure 7 shows 
that one obtains very good results for the real part of t^^^ with ct/ = 5 MeV, while there is 
somewhat less stabihty for the inversion results with uj = 20 MeV. The comparison with 
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the conventional calculation is also here satisfactory. There are only deviations of about 
1%. Again we list the relative FSI effect comparing with the Born term: -230% (/ = 0), 
+2.5% (/ = 2). The imaginary parts are somewhat more problematic. For I = one has 
the already mentioned problem with the elastic contribution combined with the fact that 
one needs Roo{E) close to threshold {Enp = 1 MeV), but with cr/ = 5 MeV one obtains 
a sufficiently good result as seen from the comparison to the result of the conventional 
calculation. Though the relative differences to the conventional calculation are rather large 
for the imaginary part of the / = 2 transition, its value is in principle correct, since it is 
more or less identical to 0. Note that it is about 200 times smaller than the already very 
small real FSI part of the / = 2 transition. 

We do not show results for the angular distribution of fi. However, from the discussion 
above it should be clear that the two different calculations lead for /l(^) to relative differ- 
ences of considerably less than 1% for kinematics (i) and (ii) and of about 1% for kinematics 
(iii). 

As mentioned in Sec. || one has to use a somewhat different method for the calculation 
of the LIT for an exclusive reaction with more than two fragments in the final state. Two 
other possibilities are discussed at the end of Sec. In both cases one has to solve different 
differential equations, e.g., (^) instead of Eq. (^). However, these new methods appear 
to be numerically more problematic. Small errors in the solution of the differential 
equation (0), might lead to a much larger error for the solution of Eq. (^), where \&i 
serves as source term on the rhs. Also for the determination of the LIT via Eq. (^) it is 
important how precise and are calculated, since one has to determine the difference 
{^^'^\V\^[{a)) - {^^^\V\^i{a)). We are able to study this question for the d{e,e'p)n 
reaction, since we can also use these alternative ways of evaluating the LIT. As a matter of 
fact both alternative methods lead in our case essentially to the same results with relative 
differences smaller than 0.01%. In Fig. 8 we show for a few selected cases these new LIT 
results compared to those obtained with Eq. (|TBp. On finds relative deviations of the order 
of 1%. There are larger differences for the kinematics with Enp = 1 and 120 MeV beyond 
a aji of 150 MeV, but they are rather unimportant, since both i?oo are very small there. In 
fact Rqo crosses zero at about 205 and 190 MeV for E^p = 1 and 120 MeV, respectively. 
Altogether one can say that one does not encounter greater numerical problems in evaluating 
the LIT with these alternative ways. Therefore also a calculation of an exclusive reaction to 
a three-body channel should lead to rather reliable results with the LIT method. 

V. CONCLUSION 

We have calculated the longitudinal response of the exclusive d{e, e'p) reaction with the 
method of the Lorentz integral transform. This method allows one to include the complete 
FSI, however, without explicit use of final state wave functions. It is the first time that the 
LIT method is applied to an exclusive reaction. In the past only inclusive processes have 
been studied with the LIT. The great success of the method raised the question whether it 
can also be successfully used in exclusive reactions. The results in this work show that one 
obtains a very precise determination of the various transition matrix elements. Differences 
to the conventional calculation are generally below 1 %. Only in the case of a very small 
FSI effect on the transition strength, i.e. a 10~^ effect compared to the corresponding Born 
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term, one can also obtain somewhat higher differences of a few percent. However, in this 
case differences could, as discussed in Sec. |^ at least partly be due to a small inexactness in 
the conventional calculation. There is only one exception, where one can expect a somewhat 
larger size of the error of the LIT result. This is the case for a transition matrix element 
in a region with transition strength from two (or more) rather narrow lying peaks. We had 
chosen such a situation with our kinematics (iii), where we have a strong elastic contribution 
at about E = —2.2 MeV and another strong peak right above breakup threshold. In such a 
situation one should try to improve the resolution of the transform L{aji, aj) by choosing a 
smaller aj. In fact our results improve significantly from cx/ = 20 MeV to o"/ = 5 MeV. 

In case of an exclusive reaction with more than two fragments in the final state one 
cannot proceed exactly in the same way as for the breakup in two fragments. In this case 
one has to use other ways for the determination of the LIT. We could show that also these 
alternative methods lead to rather precise results. Therefore, in general, we may conclude 
that the LIT method leads to reliable results not only for inclusive, but also for exclusive 
reactions. 
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FIGURES 

FIG. 1. Lorentz integral transforms Ra{aji,ai = 20 MeV) for / = (top) and I = 2 (bottom) 
at g^.m. = 5 fm^^ for various Enp as shown in the figure. 



FIG. 2. Inversion result Rii{E) with aj = 5 MeV (left) and aj = 20 MeV (right) for / = 
(top), / = 2 (middle), / = 4 (bottom); as indicated in the figure curves are shown for a number of 
inversion basis functions = 6, 16, 18, 20. 



FIG. 3. Results for real (left) and imaginary (right) parts of inoi for kinematics (i) with / = 
(top), / = 2 (middle), / = 4 (bottom) as function of the number of inversion basis functions: 
diamonds (cj/ = 5 MeV), squares (cj/ =20 MeV); also shown are the results of a conventional 
calculation with explicit final state wave function (full curves) and deviations of ±1 % from these 
results (dashed curves). 



FIG. 4. As Fig. 2 but for kinematics (ii). 



FIG. 5. As Fig. 3 but for kinematics (ii). 



FIG. 6. As Fig. 2 but for kinematics (iii) and Z = 0, 2. 



FIG. 7. As Fig. 3 but for kinematics (iii) and Z = 0, 2. 



FIG. 8. Ratio of Lorentz integral transforms i?oo(c"i?) = 5 MeV) calculated with Eq. (21) 
(dots) relative to the results of Eq. (|l|) at ql^ = 5 fm'^ for Enp = 1 MeV (top), 120 MeV 
(middle), 200 MeV (bottom). 
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